High temperature thermal conductivity of 2-leg spin-1/2 ladders 
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Based on numerical simulations, a study of the high temperature, finite frequency, thermal con- 
ductivity K,[u)) of spin-1/2 ladders is presented. The exact diagonalization and a novel Lanczos 
technique are employed. The conductivity spectra, analyzed as a function of rung coupling, point 
to a non-diverging dc— limit but to an unconventional low frequency behavior. The results are dis- 
cussed with perspective recent experiments indicating a significant magnetic contribution to the 
energy transport in quasi-one dimensional compounds. 

PACS numbers: 66.70.-|-f, 75.10.Pq, 75.40.Gb, ,75.40.Mg 



Introduction — Recent experiments Q, H, H, ^ 1^ con- 
vincingly promoted magnetic excitations as a very effi- 
cient mechanism for energy transport in quasi-one di- 
mensional materials. In particular, spin-1/2 Heisenberg 
chain and ladder compounds (undoped or hole doped) 
have been studied. In both systems, the observed highly 
anisotropic thermal conductivity (comparable in magni- 
tude to that of metallic systems) was attributed to mag- 
netic transport. 

Theoretically, it has been noticed that in the one 
dimensional (ID) spin-1/2 Heisenberg model the en- 
ergy current commutes with the Hamiltonian and thus 
the thermal conductivity is ballistic at all temperatures 
0, 0, ■ This observation falls in line with a proposal 
of unconventional transport in ID integrable systems Q . 
Besides this rigorous result, the experimental develop- 
ments motivated theoretical studies of the thermal con- 
ductivity in other (e.g. frustrated, gapped) spin chain 
and ladder Hamiltonians based on numerical simulations 
0, Hol lll| or low energy effective theories [ll[l3. These 
works focused on the thermal Drude weight Dth as a cri- 
terion of baUistic transport 0, Q, 0| and they lead to 
an active discussion whether nonintegrable systems also 
show ideal thermal conductivity. 

The key to obtaining reliable information by numerical 
simulations is the study of large enough systems that al- 
low analysis of the scaling with lattice size. In this work, 
we study the thermal conductivity within linear response 
theory, first, at high temperatures where numerical simu- 
lations on finite systems are most reliable (the character- 
istic thermal scattering length should be smaller than the 
size of the lattice) and second, by using a novel technique 
based on the Lanczos method and the microcanonical en- 
semble (MCLM [3). This method gives access to fairly 
larger systems than those studied till now by the exact 
diagonalization (ED) method. 

The main issues of actual experimental and theoreti- 
cal interest that we address are: (i) whether the thermal 
transport is ballistic or diffusive, (ii) how the interchain 
coupling affects the ideal thermal transport of the decou- 
pled (integrable) Heisenberg spin-1/2 chains, (iii) what is 



the order of magnitude of the magnetic contribution to 
the thermal conductivity as a function of exchange cou- 
plings and in particular how it compares with experimen- 
tal values in ladder compounds. 

Hamiltonian and Method — The 2-leg ladder Hamilto- 
nian is given by the g = component of Hq, the Fourier 
transform of the local energy density, that we define as, 

Hq ^ e*«('+^/'nSi.i+i ■ Si., + S2,,+i • S2,i) 



J 



(1) 



S/ are spin-1/2 operators at site I and in the following, 
we consider systems with periodic boundary conditions, 
we take J = 1 as the unit of energy (7i = fcs = 1) and 
unit lattice constants. 

To study the thermal conductivity we define the energy 
current operator using the continuity equation for Hq 
[itIIT^ . Extracting the long- wavelength transport limit, 
dHq/dt ^ —iqj^ for q ~^ 0, we obtain. 



(2) 



j'' = r^i i^u-i ■ (Si,i X Si,,+i) 

1=1, L 

+ S2J-I • (S2J X S2,,+l)} 

+ lyiSi,;-! • (Si,z X 82,0 -I- S2,z-i • (82,, X Si,0 

+ S2.1 ■ (Si,, X Si,+i) + Si, • (S2,, X S2,,+l)}. 

Within linear response theory [l^ , the real part of the 
thermal conductivity at frequency uj is given by. 



k{uj) = 2'KDthS{uj) + Kregiuj), 

with the regular part. 



(3) 



Kreg{uJ > 0) = — -tanh(-— )Si / dt 
LoL 2 Jo 



'{{f{t).f}), 
(4) 



{(3 ~ 1/T, z = Lu + irf) and the thermal Drude weight, 

DtH^^y.Pn\{n\f\n)\''. (5) 
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I n) are the eigenstates and p„ the Boltzmann weights. 
For /3 — > we can deduce the sum-rule, 



(6) 



that suggests the analysis of 2TTDth //in order to estimate 
the contribution of ballistic transport. 

In the following, we employ the ED method to study 
Dth and k{uj) for systems up to L = 9 rungs. For larger 
systems, up to L = 14, we use the MCLM method. In 
both cases we implement the translational symmetry that 
provides independent A:— Hamiltonian subspaces. In the 
MCLM method 0|, we replace the thermal average in 
Q by the expectation value over a single state | A) of 
energy A that equals the canonical ensemble value of the 
energy at the desired temperature |23|, A = (H). In 
most cases, we present results for /3 — > where A ^ 
due to the symmetric spectrum of the Hamiltonian. 
Otherwise, we determine A by extrapolation of thermal 
energies evaluated using ED results. 

The state | A) is constructed by employing a 1st Lanc- 
zos procedure of about 1000 steps using as "effective 
Hamiltonian" the operator K = {H — A)^; in practice, 
for large systems with dense spectra, the procedure can- 
not fully converge. | A), the ground state of K in the 
constructed Lanczos subspace, is characterized by a dis- 
tribution over the eigenstates | n) of variance in energy 
of 0(0.01) that imposes a maximum uj resolution to the 
spectra. Then, starting from | A), a 2nd Lanczos 
procedure of about 4000 steps provides k{uj) using the 
continued fraction technique [21| (typically r/ 0.01). 
Thus we effectively evaluate k(w) by. 



.(.) ^ ^tanh(^)(-c,)(^,|,..__l__ 



1 



z+iH-\) 



f I A)). 



(7) 



Note that, as the Lanczos procedures do not fully con- 
verge, an eventual 5(ti;)-peak contribution appears as a 
broadened weight at very low frequencies |1(t| . 

Thermal conductivity — The first step in characteriz- 
ing k{u!) is the evaluation of the Drude weight at finite T. 
If Dth is finite then the transport is ballistic; if it vanishes, 
the transport is normal provided the Kdc = ^(u; 0) 
limit exists. For a finite system Dth is always nonzero, 
so it is crucial to examine its scaling as a function of size. 
Exact results obtained by the ED method are shown in 
Figure nj In the high temperature limit Dthl0^ rapidly 
decreases, seemingly exponentially fast, with system size 
(for each series of even or odd V). It already represents 
only a couple of percent of the sum-rule for L = 9 rungs 
(inset). Of course, the size of the studied lattices is rather 
limited but a vanishing Drude weight in the /? — > limit 
is corroborated by the k(uS) spectra shown below. 
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FIG. 1: Thermal Drude weight for Jx = 1 as a function of 
temperature for systems with L rungs (ED evaluation). Inset: 
normalized Dth for /3 — > as a function of system size. 




FIG. 2: Normalized thermal conductivity for different L's by 
the MCLM method (curves displaced for clarity). ED evalua- 



i I A) tion for L = 8: thin line - k{u!)/I, thin dashed line - 



,(^)/J. 



At lower temperatures we find a non-monotonic scaling 
for even-rung systems. Note however that the change of 
scaling behavior of Dth^ from decreasing with increasing 
L at high-T to the opposite at low-T, shifts to lower 
temperature as we consider larger systems (see also [ll| ) . 
We can then argue that over the whole temperature range 
the Drude weight scales to zero as L ^ oo. It would be 
unexpected if there is a transition to a finite Dth below 
some critical temperature. 

Next, we discuss the thermal conductivity as a func- 
tion of UJ and L in an attempt to determine whether it 
represents normal transport (diffusive behavior exempli- 
fied by a Lorentzian form), unconventional (e.g. power 
law at low frequencies) or some other behavior that still 
implies a finite dc— conductivity. Furthermore, we study 
the dependence on J± in order to find out the effect of in- 
terchain coupling on the ballistic transport of decoupled 
chains. For Jj^ = only the 2'KS(uj)Dth contribution 
exists as [j^ , H] — 0. 
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In Figure|5]we show a series of spectra for different L's, 
J_L = 1 and (3^0. To illustrate the applicability of the 
MCLM method, we compare a spectrum from a canoni- 
cal ensemble ED study for L = 8 (weights smoothed with 
Lorentzians) with one obtained by the MCLM averaged 
over all /c— symmetry subspaces (in both cases 77 = 0.05). 
We find that the agreement is fair and the broadened 
Drude (5 (w)— peak is reproduced. We also show the regu- 
lar part from the ED evaluation indicating that the pres- 
ence of the Drude weight, that is fairly large for L = 8, 
is limited to w < 0.3. 

Next, in order to examine the finite size scaling of k(w), 
we present spectra for larger systems, L = 10, 12, 14 
(77 = 0.01). We notice that the statistical fluctuations 
drastically decrease with increasing L so that for L — 14 
it is sufficient to consider only one /c— subspace (fc = is 
shown, the curves being practically indistinguishable for 
the other fc— subspaces). The reason is that the dimen- 
sion of the Hilbert space increases by a factor ^10 when 
increasing L by 2 rungs. For L = 14 it is 0(3 x 10^) 
states in each fc— subspace and so it represents a fairly 
dense spectrum at high energies (temperatures). Thus, 
in the following, we show spectra constructed with only 
the fc = subspace for L = 14 rungs. 

Finally, we note that the spectra rapidly collapse to the 
same curve. Noticeable deviations for L > 10 are limited 
to the w ~ — 0.3 range. The data apparently converge 
to a K(tt') curve that: (i) it is characterized by a finite 
Kdc = k{u! — > 0) value, (ii) it approaches the uj —>■ limit 
with a finite slope, (iii) it shows a change of curvature 
at UJ ^ 0.5 for all L, (iv) its overall oj dependence does 
not correspond to a Lorentzian or power law behavior; 
a minimal description of k{uj) is obtained by an e"!'^''^ 
form [13. 

Of course, the deduced w — > behavior is rather ten- 
tative considering the limited size of the systems that we 
are able to study. A finite L imposes a cutoff on the 
lowest frequency behavior that can be reliably extracted. 
We cannot exclude the emergence of a more conventional 
form - with vanishing zero frequency slope - in the range 
u! ^ — 0.2. The apparent finite slope might also be due 
to the mixing, in the very low frequency spectrum, of a 
remnant (5(cj)-Drude contribution that we expect to diss- 
apear as L — > 00 (see Figure^. In any case, a prominent 
Drude peak is not observed. 

Next, in Figure |21 we show a series of spectra for 
L = 14 as a function of interchain coupling J± in the 
/3 — > limit. For J± < 1, first, the c?c— conductivity (in- 
set) scales as Kdc ^ ^/J'±: second, k{uj) correctly tends to 
a (5 (w)— peak as J_l — > signafing the ballistic transport 
of decoupled chains, third, the frequency range of finite 
weight extends up to iv ^ 4:J±. For Ji^ > 1, a second 
weak peak appears a.t uj 2J± (corresponding to the en- 
ergy of singlet-triplet rung excitations) and Kdc ~ J± as 
now J_L becomes now the dominant energy scale. Notice 
that for J± < 1, the overall behavior of k{uj) is what one 
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FIG. 3: (a) Thermal conductivity for L = 14 rungs as a 
function of J± for /3 ^ 0. Inset: dc— limit. 




FIG. 4: Thermal conductivity for L = 14 at T = 1, T = 
2(/3 — 0.5) and for comparison at /3 ^ 0. 

would qualitatively expect by taking into account within 
a perturbative scheme the effect of interchain coupling 
on the ballistic transport of a non-interacting system. In 
this case however, the "free" limit is the spin- 1/2 Heisen- 
berg Hamiltonian which implies that the effect of a per- 
turbation on the ideal transport of an integrable system 
is described within a simple picture. 

In Figure 01 we show the frequency dependence of the 
thermal conductivity at lower temperatures. An over- 
all similar behavior as for /3 ^ is found except that 
the rounding at a; 0.5 disappears for T = 1. The 
data are qualitatively described by two separate expo- 
nential forms, one in the frequency range w ^ — 0.7 
and another one above. No excessive enhancement of 
the (ic— conductivity is observed. The sum-rule jnj is 
now satisfied to about 95 percent and the L dependence 
of the conductivity spectra remains within the statistical 
noise of the data. It is rather difficult to study even lower 
temperatures as the energy spectrum is very sparse at low 
energies and thus the fluctuations signiflcantly increase. 

From the extracted dependence of Kdc on J± we can 
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estimate an experimentally relevant 



magnetic 



contribution to the thermal conductivity at high temper- 
atures. Re- inserting units in the basic expression Q), 



^exp 



is given by, 



^exp 



a, 6, c are lattice constants (c along the ladder axis) and 
energies are in degrees K. For characteristic lattice con- 
stants of O(lOl) and J ^ O(IOOOK), we obtain a ther- 
mal conductivity 0(2 W/mK) that scales as (/3J)^ and 
( J/ J^)^ for T 3> J. If we assume an exponential increase 
of the conductivity (due to the freezing of Umklapp pro- 
cesses) with characteristic energy J below a temperature 
of 0( J), we obtain a room temperature thermal conduc- 
tivity of order 10 - 100 W/mK that is consistent with 
experiment. We should keep in mind in this estimation 
that it is notoriously difficult to establish the temper- 
ature dependence of the thermal conductivity for T < J. 

We can then argue, that the experimentally observed 
high values of the thermal conductivity 0, Q in ladder 
compounds arc due to the large exchange coupling J and 
that they are limited by spin-spin as well as spin-phonon 
scattering. This is in contrast to quasi-lD materials de- 
scribed by the Heiscnbcrg spin- 1/2 Hamiltonian [3, 
where it is unambiguous that the spin-spin scattering is 
absent. 

Discussion - The presented numerical study indicates 
that, at least in the high temperature limit, the ther- 
mal conductivity in 2-leg spin- 1/2 ladders is character- 
ized by a finite c?c— value and thus non ballistic trans- 
port. The low frequency spectra do not seem consistent 
with a Lorentzian form. At lower temperatures there is 
a change in behavior but not a dramatic enhancement 
of the dc— conductivity or the appearance of a promi- 
nent Drude peak. Presently, with the available numerical 
methods, it is rather challenging to explore an eventual 
crossover to a quasi-ballistic regime at low temperatures. 

Furthermore, this analysis exemplifies the effect of a 
non-integrable interaction (interchain coupling) on the 
ballistic transport of an integrable system. The ob- 
tained perturbative result indicates that the presence in 
a Hamiltonian of an integrable case with diverging con- 
ductivity is signaled over a finite range in interaction pa- 
rameter space. 

Regarding experimental realizations, this study is more 
relevant to ladder compounds with J smaller than the 
room temperature but it shows that materials with high 
thermal conductivities can be obtained by increasing the 
coupling J and/or reducing the interchain coupling J±. 
It also shows that it would be interesting to explore the 
unusual frequency dependence, for instance, by light scat- 
tering. As to the observability of the exceptional trans- 



port in integrable systems, we expect the effect to be 
most remarkable at high temperatures in small exchange 
coupling J one dimensional spin-1/2 Heisenberg systems. 
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